Electronic structure of CdTe nanocrystals: A tight-binding study 
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We present a symmetry-based calculation of the electronic structure of a compound semicon- 

^ . ductor quantum dot (QD) in the sp'^s* tight-binding model including the spin-orbit interaction. 

^ ' The Hamiltonian matrix is diagonalized exactly for CdTe QD sizes up to 60 A. The surface dan- 

[L , gling bonds are passivated by hydrogen through a careful analysis of the density of states and wave 

functions. The calculated size dependence of the energy gap shows a reasonable agreement with 

the available experimental data. Our symmetry analysis indicates that, in contrast with a reported 

prediction of the three-band effective-mass model, the fundamental interband transition remains 

dip ole- allowed in CdTe nanocrystals. 

C^ 

^ I PACS Numbers: 73.20.Dx, 71.24.-|-q, 71.15.Fv, 78.40. Fy 
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I. INTRODUCTION 



Semiconductor nanocrystallites, also called suantum 
dots (QD's), have been studied in recent years.EI In par- 
' O ticular, high-quality crystallites of II- VI compounds have 
^ been fabpinated and their optical properties investigated 

7^ I in detail.aEl The confinement-induced blueshift of the fun- 
damental gap and the discretization of the energy sp' 
trum have been observed. Early theoretical studiei 
were based on the effective mass approximation (EM A), 
which was progressively refined byrjtaking into account 
the degeneracy of the valence band.ffQ Generally speak- 
,_^ ' ing, the EMA overestimates the ccpfinement energy in 
psj > small QD's. Lippens and Lannootd reported a tight- 
^^ ' binding (TB) calculation of the QD energy gap based 
On on the ernpirical sp^s* model, previously developed by 
^^ Vogl et aZ.til for bulk semiconductors. They used the re- 
cursion method and obtained a better agreement with cj 
perimental data than the EMA calculations. Nair et al 
presented a diffcrent-lTB approach based on an effective- 
bond-orbital modehlla The pseudopotential-oicthod has 
^ I also been succesfciUy applied to some QD's.E£rEZl 
O . Ren and Dowlla developed an algorithm for exact diag- 

O ■ onalization of the TB Hamiltonian for a hydrogenated Si 
J> \ cluster in the sp^s* model by using tiie tetrahedral (T^) 
symmetry. More recently, Albe et aZJla reported TB cal- 
culations for ZnS and CdSe QD's, also based on an exact 



J3 ■ diagonalization of the cluster Hamiltonian, but including 

" " ' the spin-orbit interaction. However, the "brute force" 

method, suitable for studying the QD shape effects they 

focused on, is rather limited to small sizes and does not 

allow any symmetry analysis of the QD eigenstates. 

In this paper we report a symmetry-based TB study 
of the electronic structure of a zincblende semiconductor 
crystallite of roughly spherical shape in the sp^s* model 
including the spin-orbit interaction. Thus, our approach 
generalizes that of Ref. |l^ to binary compounds and a 



finite spin-orbit interaction. Here wei-apjaLv it to CdTe, 
because substantial experimental datal3t3'£d are available 
and a thorough three-band EMA studyQ has been re- 
ported, allowing a detailed comparison with our results. 
The paper is organized as follows. We first present an 
outline of the model in Section II, including the method of 
passivation of the surface states through hydrogen bond- 
ing. We present and discuss some results for CdTe QD's 
in Section III. The size dependence of the energy and 
the symmetry classification of a few conduction and va- 
lence band edge states are shown. The variation of the 
energy gap is compared with the available experimental 
data and the results of other pubhshed calculations. Fi- 
nally, an analysis of the orbital symmetry of the highest 
valence band and the lowest conduction band states is 
carried out in order to check the electric dipole selection 
rule for the fundamental transition. The paper closes 
with some concluding remarks. 



II. THE MODEL 

The smnpiiidependent part of the TB Hamiltonian is 
given by,ElO 

Ho= ^ \Rb,i,fT)Ei^t(Rb,i,(j\+ ^ |Ra,i,cr)K:j(Rc,j, cr| -l-H.c. 

(1) 

where b = a (anion), c (cation). R(, are the atomic po- 
sition vectors, i and j denote the orthonormal atomic 
orbitals s, Px, Py, Pz and s*, which depend on b. a is the 
z component of spin (| and j). Note that s* represents an 
an excited s state introduced phenomenologically in order 
to obtain a correct description of the conduction band in 
bulk semiconductors. The interatomic matrix elements 
Vij are restricted to the nearest neighbors as indicated 
by < Ra,Rc > in the summation index. 13 indepen- 
dent parameters then characterize Hq for a compound 
semiconductor, which are chosen to fit the known band 
structure— For example, the parameters for CdTe are 
given byEl Es,a = ~8.891, Ep^a = 0.915, Es,c = -0.589, 
Ep,c = 4.315, Vs,s = -4.779, 14,, = 2.355, Vx,y = 4.124, 
Vs,p = 1.739, Vp.s = -4.767, Es..a = 7.0, S,.,e = 7.5, 
Fs*,p = 1-949, Vp^s' = -2.649 eV. The spin-orbit cou- 
pling part of the Hamiltonian mixes the spin- up and spin- 
down p orbitals on the same atom: 

Hso= J2 \Rb,i,T)2XbU-Sb{Rb,j,'y'\- (2) 

It introduces 2 additional parameters. \a = 0.367 and 
Ac = 0.013 eV for CdTe.ta When Hso is included a more 
convenient atomic basis is given by the spin-orbit coupled 
orbitals, which are basically the total angular momentun 
eigenstates |Rb, j, jz) (j = 3/2, 1/2) that diagonalize the 
spin-orbit interaction. We, however, write them in terms 



of the basis functions |Rb,u^„) of the irreducible repre- 
sentations Ffc (fc = 6, 7, 8) of the tetrahedral point group 
Td with respect to the site Rf,. Note that we have two 
sets of Fg orbitals at a given site arising from the s and 
s* states. Using the phase convention and couphng coef- 
ficients of Ref. |l5|, the first rows are 

l"'-3/2(p)) - -i-,\P^)\i) + ^,\Py)\i) + ^IP^)IT)- 

Let us now construct a cluster of roughly spherical 
shape starting from, say, a cation at the origin by suc- 
cessively adding nearest-neighbor atoms through tetra- 
hedral bonding. The dangling bonds emanating from 
the atoms generated in the last step will be passivated 
by placing a hydrogen s orbital at each empty nearest- 
neighbor site. This "hydrogenated" crystallite has an 
overall tetrahedral symmetry. We can, therefore, reduce 
the Hamiltonian to a block diagonal form by rewriting 
it in a symmetrized basis corresponding to the double- 
valued representations F^ (fc = 6, 7, 8) of T^. 

Following the procedure developed in Ref. |lj, we 
first construct the symmetrized site functions (j)\, cor- 
responding to the single- valued representations Ti{i = 
1, 2, 3, 4, 5) of Td- All atoms at a given distance Rb from 
the origin constitute a shell. The sites {Rb} on a given 
shell are grouped into (generally more than one) sym- 
metry subshells, each containing all sites which trans- 
form into one another under the 24 symmetry opera- 
tions of Td with respect to the QD center. Each subshell 
{R} is then spanned by the functions 0^ which repre- 
sent symmetrized linear combinations of sites, assuming 
the values 0^(R) at R. These values are deduced by 
using the projection operators. Note that, except for 
the one-dimensional representations Fi and F2, an irre- 
ducible representation generally occurs more than once 
in a given subshell. The site functions are then coupled 
with the localized atomic orbitals and the subshell basis 
functions of total symmetry Ffc(fc — 6,7,8) obtained as 
linear combinations of ^j^(?!)5j(R)|R, u^i), by using the 
coupling coefficients for Tk in F^ x Tj. 

The hydrogen shells are treated in a similar way. Here 
we have only one s-orbital per atom so that the local ba- 
sis is restricted to the Fg symmetry: IRfjiMm}. Finally, 
in the symmetrized basis, the total Hamiltonian becomes 
block diagonal, each block corresponding to a given sym- 
metry Ti{i = 6, 7, 8). This implies a substantial reduction 
of the size of the matrix to diagonalize. For example, for a 
cluster containing 3109 semiconductor atoms plus 852 H, 
instead of a matrix of size 32794 = 10 x 3109 -I- 2 x 852 we 
get two equivalent Fg matrices of size 2740, two equiv- 
alent F7 matrices of size 2739 and four equivalent Fg 
matrices of size 5459. Thus, the matrix size has been 
reduced by a factor of 6. 

While the TB parameters enumerated above are as- 
sumed to be the same as those for the bulk semiconduc- 
tor, we need to precise the parameters concerning the hy- 



drogen atoms. The H energy level is obtained through the 
same scaling prescription as that for the cation and anion 
s levels.li3 In the case of CdTe, we obtain Eg^n — —5 eV. 
The hopping matrix elements between the anion or cation 
and H are assumed to follow the Harrison scaling rule: 
Vf,_H = {da-c/ dh-iii)^Vac, in terms of the bond lengths. 
dh_H are usuallyll3 taken from molecular data. In the ab- 
sence of such data on Cd-H or Te-H bonds, we estimate 
the bond lengths by summing the corresponding covalent 
radii. This procedure typically yields the experimental 
bond lengths to within a few hundredths of angstrom. 
We thus obtain dcd-H = 1-71 and dxc-H — 1-67 A. The 
above parametrization scheme is adequate for passivatiBS 
the surface dangling bonds in Si and Ge nanocrystals.tS 
However, in the case of CdTe, as shown in the next sec- 
tion, it is necessary to readjust (reduce) the bond lengths 
in order to completely eliminate the surface states from 
the electronic spectrum in the neighborhood of the for- 
bidden gap. 



III. RESULTS AND DISCUSSION 

Let us first discuss the saturation of dangling bonds, 
as we are basically dealing with surface-passivated QD 
systems which show a systematic increase of the opti- 
cal energy gap with decreasing size. The electronic den- 
sity of states (DOS) as well as its hydrogen projected 
part, as calculated from the exact eigenstates of the TB 
Hamiltonian, are partly shown in Fig. 1 for a CdTe QD 
of diameter D — 58.67 A containing N — 3109 semi- 
conductor atoms. Note that D = a(3iV/47r)i/3. The 
uppermost spectrum (a) corresponds to the CdTe clus- 
ter without H. Notice the large number of states in the 
bulk band gap region between and 1.6 eV, which arise 
from the surface dangling bonds. Next we show the re- 
sults (b) for the hydrogenated cluster (852 H atoms) with 
the standard bond length parameters: dcd-H = 1-71 and 
dxe-H = 1-67 A. Most of the surface states have disap- 
peared, but not all. In fact, as seen from the hydrogenic 
partial DOS superposed on the total DOS, the lowest un- 
occupied (conduction) states are still surface states. We 
have directly verified this through the spatial distribution 
of the wave function. This is in contrast with the cases of 
Si or Ge where a complete saturation of dangling bonds 
is achieved by using the molecular bond lengths for the 
couplings with H atoms. As shown in Fig. 1 (c), in order 
to eliminate surface states from the relevant part of the 
electronic spectrum in our CdTe QD, we need to increase 
the couplings with H. By choosing bond lenghts slightly 
shorter than before, dcd-H = 1-58 and dTc-H = 1-54 A, 
the hydrogenic DOS in the conduction band vanishes up 
to ~ 2 eV. We have checked that the same set of bond 
length parameters assures surface passivation in all the 
QD's considered. 

Now we present the energies and wave functions of 
CdTe quantum dots of several sizes (see Table B). We 



have chosen cation- and anion-terminated crystaUites al- 
ternatingly in order to show the typical oscillating behajv-, 
ior of energy levels in the clusters of binary compounds.t3 
For instance, two crystallites of similar diameters, 12.35 
and 13.15 A, and terminated by cation and anion shells, 
respectively, present the lowest unoccupied state (LUS) 
at 2.11 and 2.42 eV. This is in contrast with the 
monotonous decrease observed in Si or Ge crystallites. 
The valence states also show such oscillations but with 
a smaller amplitude. In general, the LUS energy of a 
cation-terminated crystallite is higher than that of an 
anion-terminated crystallite of neighboring size, because 
of the dominating contribution from the cation s orbital 
to the conduction band. In Figures ^ and ^ we plot a few 
energy levels below and above the bandgap, respectively, 
against the QD size. The symmetry classification of the 
levels is also indicated. Note that for all sizes the LUS 
(conduction band) is of Fg symmetry and the highest oc- 
cupied state (HOS) of the valence band belongs to Fg, as 
in the bulk semiconductor. For comparison, we include 
the respective energies calculated in the EMA (Ref. ||) 
which are shown as the solid curves. As usual, the EMA 
confinement energies are much larger than the TB ones 
in small-size QD's. However, the HOS (Fig. 2) shows 
a better convergence with increasing size, which reminds 
us that the TB model yields a better description of the 
valence band in the bulk semiconductor. 

The calculated QD bandgap Eg, the energy difference 
between the LUS and HOS, is plotted as a function of 
size in Fig. ^. We also present a compilation of the avail- 
able experimental values of the optical bandgap. The re- 
sults of previous EMA (Ref. |) and TB (Ref. ]l^) calcula- 
tions are also shown. In contrast with these calculations, 
however, we have not corrected Eg for the electron-hole 
Coulomb interaction. Thus, tis.e closeness of our results to 
those of Lippens and Lannoot2l shown as the dashed curve 
is rather misleading. Their uncorrected gap would be sys- 
tematically larger than ours; this is probably related to 
the neglect of spin-orbit coupling in their calculation. It 
can be seen that the size dependence of our uncorrected 
Eg shows a good qualitative agreement with experiment. 
The calculated values are somewhat smaller, but they are 
simply based on the parameters of Ref. |l2|for bulk CdTe, 
which we prefer not to modify in any arbitrary manner. 
There is also some uncertainty related to the parameters 
concerning the H atoms. As for the Coulomb correction, 
it is rather difficult to evaluate consistently within the 
TB model. But we expect it to be significantly smaller 
than that estimated from the EMA used in the calcula- 
tions cited above, because the TB carrier wave functions 
are spatially more extended (see Fig. 5). 

Figure 5 shows the shell- wise radial distribution of the 
carrier probability density for the LUS (electron) and 
HOS (hole) states in a QD of diameter 58.67 A. The 
shell probability density for an eigenstate |^) is directly 
obtained from the diagonalization of the TB Hamilto- 
nian, by summing the local density P(Rfc) = ^^ ^ | < 



RbjijCrl^f > p over the sites Rb on the given shell. We 
see that, in contrast with the smooth EMA envelope func- 
tions, the TB radial distributions are oscillatory. It is in- 
teresting to note that the LUS (HOS) is indeed prcpon- 
derent at the cation (anion) shells. In the present case 
the cation (anion) shells account for 79% (76%) of the 
LUS (HOS). Another reason for the oscillatory behavior 
is the discontinuous variation of the number of atoms on 
a shell. All the oscillations are, of course, smoothed over 
in the EMA. It is not easy to make a direct comparison. 
However, it can be verified that our TB charge distri- 
butions have a greater radial extension than the EMA 
envelope functions. 

Finally, in order to reexamine the electric-dipole selec- 
tion rules for the fundamental interband transition, we 
present an analysis of the orbital symmetry of the full 
wave functions in the TB model. We proceed as follows. 
By setting the spin-orbit interaction -ffso = 0, we first di- 
agonalize Hq in the spin-degenerate atomic orbital basis. 
The resulting eigenvalues are then classified according to 
the single- valued representations F^ (j = 1 — 5). As the 
second step, we diagonalize the same Hamiltonian in the 
spin-orbit coupled basis: the same eigenvalues are now 
classified in terms of the double-valued representations 
Fi {i — 6 — 8), which arise from the direct products of 
Fg with Ti (i = 1 — 5). By consulting the multiplica- 
tion table of the group Td, we see that there are three 
different types of Fg corresponding to the three orbital 
symmetries F3, F4 and F5. Thus, the eigenstates Fg of 
Ho can be identified as Fg(F3), Fg(F4) or Fg(F5). Finally, 
a given eigenstate Fg of the full Hamiltonian, for example 
the HOS, can be written as a linear combination of all 
the 'unperturbed' states Fg(Fi) (i = 3 — 5). By calculat- 
ing the projections and summing the probabilities over 
all states arising from a given orbital symmetry F^, we 
can deduce the fractional contribution of that symmetry 
in the HOS. Similarly, an eigenstate Fg such as the LUS 
can be analyzed in terms of the orbital symmetries Fi 
and F4. 

The numerical results for the CdTe QD's are as fol- 
lows. The LUS is found to be almost pure Fi for all sizes: 
99.9% for D = 12.35 A and 99.7% for D = 58.67 A. The 
results for the HOS (Fg) are more interesting; they are 
shown in Table O. Note that our numerical calculation 
was limited to the contributions from the nine topmost 
'unperturbed' valence states for each of the single- valued 
representations, because of the huge size of the data files 
involved. This is, however, adequate for identifying the 
majority component in all cases. We find that the F3 con- 
tribution is always small. On the other hand, as the QD 
size increases the dominant orbital symmetry of the HOS 
changes from F4 to F5. Thus, for D > 17.8 A the funda- 
mental interband transition is certainly dipole-allowed. 
In Fig. 2 we notice that the second highest valence state 
is also of Fg symmetry and lies very close to the HOS. 
We, therefore, analyzed its orbital symmetry as well; the 
results are presented inside parentheses in Table ^ It 
is interesting to observe that its dominant orbital sym- 



metry also changes with increasing size, but from Fs to 
r4, in a manner opposite to the HOS. Anyway, the QD 
size where this second level approaches the HOS is large 
enough so that no further 'symmetry cross-over' is ex- 
pected, because the HOS in the bulk is indeed Fs-like. 
To sum up, the fundamental transition is almost dipole- 
forbidden for D < 13.2 A, but certainly dipole-allowed 
for D > 17.8 A. Clearly, our results contradict the predic- 
tion of Lefevre et aln. They found that in the three-band 
EMA the HOS (hole ground state) in a CdTe QD, regard- 
less of the size, has a p-type envelope function, making 
it dipolc-forbiddcn for optical transitions with the s-type 
LUS. Let us point out that although our analysis is inad- 
equate for calculating the actual transition probabilities, 
it does allow an unambiguous identification of the sym- 
metry forbidden transitions. 



IV. CONCLUDING REMARKS 

We have presented a TB calculation of the electronic 
structure of a roughly spherical semiconductor QD of 
tetrahedral sypimetry It is based on the sp'^s* empir- 
ical TB modeO that accounts for the band structure of 
bulk semiconductors. The use of group theory allows not 
only a substantial reduction in size of the Hamiltonian 
matrix to diagonalize, but also a symmetry analysis of 
the eigenstates. Our treatment represents a_generaliza- 
tion of the previous work by Ren and DowtS to binary 
compounds and a finite spin-orbit interaction. Here we 
have applied the model to the case of CdTe. Through a 
careful analysis of the DOS spectrum we have explicitly 
passivated the surface dangling-bond states. The calcu- 
lated one-electron bandgap as a function of size shows 
a reasonable agreement with the available experimental 
values of the optical bandgap. We find that, regardless of 
the QD size, the valence band HOS and the conduction 
band LUS belong to Fg and Fg symmetries, respectively, 
as in the bulk semiconductor. However, an analysis of the 
orbital symmetry reveals that, while the LUS is almost 
pure Fi-like for all sizes, the dominant orbital symmetry 
of the HOS changes from F4 for D < 13.2 A to F5 for 
D > 17.8 A. Thus, except for very small QD's, the elec- 
tric dipole transitions remain allowed between the HOS 
and the LUS, in contradiction with the three-band EMA 
result reported previously. 
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TABLE I. Some important features of the QD's analyzed. 
The first column gives the diameter. The second and third 
columns give the numbers of semiconductor and hydrogen 
atoms, respectively. The fourth column specifies the atom 
type of the terminating shell. The last two columns show the 
calculated HOS and LUS energies in eV. 
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TABLE IL Size dependence of the percentage contribu- 
tions of different orbital symmetries in the HOS. Inside paren- 
theses we show the corresponding values for the nearest state 
in energy (also Fg ), which is almost degenerate for D > 30 

A. 



FIG. 1. Total density of states (dashed curve) and its 
hydrogenic part (solid curve) for D — 58.67 A: (a) No passi- 
vation, (b) H passivation with the bond lengths rfxc-H = 1-67 
and dcd-H = 1-71 A, (c) H passivation with the bond lengths 
dTe-H = 1.54 and dcd-H = 1-58 A 



FIG. 2. Energies and symmetries of the ten highest occu- 
pied (valence band) levels versus the QD diameter. The solid 
curve shows the three-band EMA result from Ref. M for the 
highest level. 



FIG. 3. Energies and symmetries of the ten lowest unoccu- 
pied (conduction) levels versus the QD diameter. The solid 
curve is the EMA result for the lowest level. 



FIG. 4. Size dependence of the energy gap: Results of the 
present calculation are shown as open triangles connected by 
solid lines. Experimental data for the optical gap are from 
Refs. pel (open diamonds), Bll (filled squares), and M ( open 
circles). Previous calculations in the TB model (Ref. [10|) and 
the EMA (Ref. H) are shown as the dashed and dotted curves, 
respectively. 



FIG. 5. Shell-wise radial distribution of the carrier prob- 
ability density in the HOS (a) and LUS (b) for D = 58.67 
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